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Abstract 

Magnetism at the nanoscale has been a very active research area in the past decades, because of its novel fundamental physics and 
exciting potential applications. We have recently performed an ab intio study of the structural, electronic and magnetic properties 
of all 3(i transition metal (TM) freestanding atomic chains and found that Fe and Ni nanowires have a giant magnetic anisotropy 
energy (MAE), indicating that these nanowires would have applications in high density magnetic data storages. In this paper, 
|we perform density functional calculations for the Fe, Co and Ni linear atomic chains on Cu(OOl) surface within the generalized 
gradient approximation, in order to investigate how the substrates would affect the magnetic properties of the nanowires. We find 
that Fe, Co and Ni linear chains on Cu(OOl) surface still have a stable or metastable ferromagnetic state. When spin-orbit coupling 
(SOC) is included, the spin magnetic moments remain almost unchanged, due to the weakness of SOC in 2id TM chains, whilst 
significant orbital magnetic moments appear and also are direction-dependent. Finally, we find that the MAE for Fe, and Co remains 
large, i.e., being not much affected by the presence of Cu substrate. 

Keywords: magnetocrystalline anisotropy, transition metal nanowires, spin-orbit coupling 
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1. Introduction 



In recent years, nanostructured magnetic materials have re- 
ceived much attention, due to their interesting physical proper- 
ties and potential applications. For example, finite free-standing 
gold atomic chains were first reported in 1998[lj 2J, and their 
structural properties, such as the actual length of the chain have 
been the focus of intensive experiments and theoretical studies 
since then. However, these free-standing atomic chains are un- 
stable and thus can only exist at low temperatures and on a suit- 
ably chosen substrate. Physically stable magnetic nanowires 
deposited on metallic substrates are one of the most important 
nanostructures and a variety of techniques have been used to 
prepare and study them. For example, Gambardella et al 01 
succeeded in preparing a high density of parallel atomic chains 
along steps by growing Co on a high-purity Pt (997) vicinal 
surface in a narrow temperature range of 10~20 K. The mag- 
netism of the Co wires was also investigated by the x-ray mag- 
netic circular dichroism.fl Experimentally, copper and tung- 
sten are excellent substrates for growth of Fe thin films [fslS] be- 
cause of the small lattice constant mismatch between Fe and Cu 
(3.61A) as well as W (3.61A). Thus, in this paper, we perform 
first principles calculations to study the electronic and magnetic 
properties of the Fe, Co and Ni nanowires. The spin-orbit cou- 
pling (SOC) is included in this study to determine the magnetic 
anisotropy energy (MAE). 



2. Theory and Computational Method 

In the present calculations, we use the accurate frozen-core 
full-potential projector augmented-wave (PAW) method, 
as implemented in the Vienna ab initio simulation package 
(VASP) (^J^. The calculations are based on density functional 
theory with the generalized gradient approximation (GGA)[ lOJ. 
A large plane-wave cutoff energy of 350 eV is used for all the 
systems considered. The Fe, Co, and Ni nanowires along the 
X direction on the Cu (001) surface are modeled by a nanowire 
attached to both sides of a seven-layer-thick Cu (001) slab as 
plotted in Fig. [T] The transition metal (TM) atoms on the 
nanowires are placed either on the top of surface Cu atoms 
[denoted here as the atop (A) site] or at the hollow position 
on the Cu surface [called here as hollow (H) site]. The two- 
dimensional unit cell is chosen to be of p(Ax\) structure. The 
nearest in-plane (out of plane) wire-wire distance is larger than 
lOA (11 A) which is wide enough to decouple the neighboring 
wires. The theoretical lattice constant (3.60 A) of bulk copper, 
which is in good agreement with experimental Cu lattice con- 
stant of 3.61 A, is used as the fixed in-plane lattice constant of 
the Cu slab. However, the atoms are allowed to move in the out 
of plane direction, and the structural relaxations are performed 
using the conjugate gradient method. The equilibrium struc- 
ture is obtained when all the forces acting on the atoms and the 
stress are less than 0.02 eV/A and 2.0 kBar, respectively. The 
F-centered Monkhorst-Pack scheme with a A:-mesh of 20 x 5 x 1 
in the full Brillouin zone (BZ), in conjunction with the Fermi- 
Dirac-smearing method with cr - 0.2 eV, is used to generate 
A;-points for the BZ integration. 
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Table 1: Calculated formation energy , equilibrium interlayer distance dcq, 
spin magnetic moment per magnetic atom, m,, and magnetization energy per 
magnetic atom ei"''s=£FM _£NM jjjg -jj jjyj n^jjowires on both the hollow 
(H) and atop (A) sites of the Cu (001) surface. Here superscripts FM and NM 
denote the ferromagnetic and nonmagnetic states, respectively. 





site 


EJ 






£mag 






(eV/u.c.) 


(A) 


(Mb) 


(eV) 


Fe 


H 


-2.14 


1.64 


3.07 


-0.337 




A 


-0.31 


2.31 


3.29 


-0.444 


Co 


H 


-2.15 


1.57 


1.79 


-0.130 




A 


-1.19 


2.27 


1.99 


-0.205 


Ni 


H 


-2.26 


1.55 


0.00 


0.000 




A 


-1.32 


2.26 


0.65 


-0.032 



Figure 1: (color online) The slab supercell for modeling Fe, Co and Ni 
nanowires on the hollow sites of Cu (001) surface. 

Because of its smallness, ab initio calculation of the MAE 
is computationally very demanding and needs to be carefully 
carried out.[ll] Here we use the total energy difference ap- 
proach rather than the widely used force theorem to determine 
the MAE, i.e., the MAE is calculated as the difference in the full 
self-consistent total energies for the two different magnetization 
directions concerned. The total energy convergence criterion is 
10"^ eV/atom. The same A:-point mesh is used for the density of 
states calculations. The MAEs calculated with a denser 32x6x1 
A;-point mesh hardly differ from that obtained with the 20x5x1 
^-point mesh (within 0.02 meV). 

3. Results and Discussion 

The chain formation energy is defined as E^ - 
j{E, - nE^""' - mEj^^") where E, is the total energy of the sys- 
tem, is the total energy of the bulk Cu, and Ej™" is the 
total energy of the freestanding transition metal nanowires. And 
n and m are the numbers of the Cu and TM atoms in the sys- 
tem, respectively. The calculated formation energy is listed in 
Table [1] Clearly, the formation energy is lower when the TM 
atoms are placed on the hollow sites, as might be expected. In- 
terestingly, among the Fe, Co, and Ni nanowires, Ni nanowire 
has the lowest formation energy on both hollow and atop sites. 
The ideal interlayer distance for Cu substrate is 1.8 A, and 
from Table [1] the interlayer distance between the Fe, Co, and 
Ni nanowires and the Cu substrate for hollow (atop) site is 1.64 
(2.31), 1.57 (2.27), and 1.55 (2.26) A, respectively. At hollow 
(atop) site, the cooper substrate seems to push (pull) the TM 
nanowires and the interlayer distance of Ni is smallest (see Ta- 
ble [TJ. Table [1] also shows that when the interlayer distance 
changes, the spin magnetic moment changes as well. In gen- 
eral, for all considered Fe, Co and Ni nanowires, an increase 
in interlayer distance will result in an increase of spin moment 
for both sites. The spin moment is larger in the atop site due 
to the smaller overlap between the TM nanowires and the Cu 
substrate. At equilibrium, the spin moments of Fe, Co and 



Ni on hollow (atop) sites are 3.07 (3.29), 1.79 (1.99) and 0.0 
(0.65) ju^/atom. In our previous studyll4] of freestanding 3t/ 
TM nanowires, the spin moments at the same bond length (2.55 
A) of Fe, Co and Ni are 3.30, 2.30 and 1.14 //g/atom, respec- 
tively. Clearly, placing the 3d TM nanowires on the hollow 
sites significantly reduces or even quenches the spin moments 
on the nanowires, whilst the spin moments are much less af- 
fected when the nanowires are deposited on the atop sites. 

The relativistic SOC is essential for the orbital magnetiza- 
tion and magnetocrystalline anisotropy in sohd, though it may 
be weak in the id TM systems. Therefore, we include the SOC 
in our further self-consistent calculations in order to study the 
magnetic anisotropy and also orbital magnetization of the Fe, 
Co and Ni nanowires on Cu (001), and the results are summa- 
rized in Table |2] First, when the SOC is taken into account, the 
spin magnetic moments for the TM nanowires at hollow (atop) 
site are 3.07 (3.28), 1.78 (1.99), and 0.02 (0.64) yUg/mag. atom 
which are almost identical to the corresponding one obtained 
without SOC. This is due to the weakness of SOC in the 3d 
transition metals. However, including the SOC gives rise to a 
significant orbital magnetic moment in TM nanowires and, im- 
portantly, allow us to determine the easy magnetization axis in 
3d nanowires. For the magnetization lies along the chain direc- 
tion, the calculated orbital moments for Fe, Co and Ni are 0. 1 0, 
0.27 and 0.02 /ig/mag. atom (see Table|2l), respectively. In our 
previous study of 3d TM nanowires, we found that the magne- 
tization has a strong directional dependence [ 14] and the orbital 
moment is larger when the magnetization lies along the chain 
direction. In this study, the orbital magnetic moments for Co 
and Ni still have a directional dependence but become smaller 
compared with freestanding nanowires, as expected, and the 
orbital magnetization is also larger when the magnetization is 
along the chain direction, except Co and Ni on the atop site. At 
the atop site, the orbital moment for Ni and Co is larger when 
the magnetization is in-plane but perpendicular to the chain di- 
rection (i.e., along the y-axis). 

The magnetic moments of the Fe, Co and Ni nanowires can 
be understood by the calculated spin-polarized density of states 
(DOS) as displayed in Fig. |2l In Fig. [2l the Fermi levels are set 
to be zero, the DOS for the minority spin are multiplied by - 1 , 
and the ip-orbital decomposed DOS are scaled by a factor of 
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Figure 2: (Color online) Spin-polarized density of states of the Fe and Co 
nanowires on Cu (001). 



Table 2: Calculated spin magnetic moment per atom, m,, orbital magnetic mo- 
ment per magnetic atom, nio, along three different directions and magnetic 
anisotropy energy per magnetic atom (MAE) of the Fe, Co and Ni chains on 
Cu (001). The chain direction is along (100), and the (010) [(001)] direction is 
in-plane [out of plane] but peipendicular to the chain direction. The MAE £' 
is defined as e'«>-E°°' and is E'«>-E°'°. 









rrio 




MAE 






(Mb) 








(meV) 








100 


010 


001 


£' 


£2 








H-site 








Fe 


3.07 


0.10 


0.09 


0.10 


0.32 


0.25 


Co 


1.78 


0.27 


0.18 


0.17 


-1.17 


-1.16 


Ni 


0.02 


0.02 


0.00 


0.00 


-0.53 


0.01 








A-site 








Fe 


3.28 


0.12 


0.11 


0.13 


0.38 


0.29 


Co 


1.99 


0.19 


0.25 


0.12 


-0.40 


-1.51 


Ni 


0.64 


0.13 


0.27 


0.11 


-0.05 


-0.32 



20 for an easy comparison. Clearly, the li-orbitals of the Fe and 
Co nanowires on both sites are significantly localized due to the 
reduction of the coordination number whilst the sp-orbitals are 
more dispersive. The reduction in coordination number also in- 
duces the great enhancement in the spin splitting of Fe and Co 
3c/-bands. The spin-splitting of Fe and Co 3t/-bands for hollow 
(atop) site are 2.58 (2.82), and 1.51 (1.75) eV, respectively. The 
interlayer distance between the TM nanowires and Cu substrate 
are larger on the atop site which means that the overlaps be- 
tween the TM and substrates are smaller, and the spin magnetic 
moments are larger The splitting of the 3ii-band is approxi- 
mately proportional to the spin moment. 

The calculated MAEs of Fe, Co and Ni nanowires are listed 
in Tabled We define the energy differences and 
£2_£i()0_gOio jf |3oth and are negative, the magnetiza- 
tion prefers to lie along the chain direction. First, Co nanowires 
have the largest MAE on both hollow and atop sites as well 
as the largest orbital magnetic moment. Interestingly, both 
and E^ are positive in the Fe nanowire, indicating that it has 
an out-of-plane anisotropy at both sites, whilst the Co and Ni 
nanowire show an in-plane anisotropy. Recent ab initio stud- 
ies of Fe double chains on Ir (001) surface 11121 11311 show that 
the anisotropy could be either in-plane or out of plane, depend- 
ing on the inter-chain distance and the magnetic configuration 
of the chain. The MAEs of the freestanding 3d TM atomic 
chains |fl4:l were already obtained by the force theorem method 
with a fine 1 x 1 xn (n = 200) ^-point mesh, and the calculated 
MAEE' was -2.25,0.68 and -8.13 (meV/mag. atom) for the Fe, 
Co and Ni chain, respectively. For a more precise comparison, 
here we recalculate the MAE of the freestanding Fe, Co and Ni 
atomic chains by the total energy difference approach with n = 
200. The resultant MAE for the Fe, Co and Ni chains is -0.6, 
0.84, -8.91 (meV/mag. atom), respectively. Clearly, when the 
3d TM chains were deposited on Cu (001), the MAE generally 
becomes smaller due to the overlap between the TM nanowires 
and substrate (Table|2]l. However, the MAE for Co is enhanced, 
though the predicted magnetic anisotropy is in-plane. Further- 



more, the MAE for the 3d TM nanowires deposited on Cu (001) 
is still very large when compared with that of bulk Fe and Ni. 

4. Conclusions 

We have performed ab initio GGA calculations for the Fe, 
Co and Ni linear atomic chains on Cu (001) surface in order 
to examine how the substrates would affect the magnetic prop- 
erties of the nanowires. We found that Fe, Co and Ni linear 
chains on Cu (001) surface still have a stable or metastable 
FM state. When the SOC is included, the spin magnetic mo- 
ments remain almost unchanged, due to the weakness of SOC 
in 3d TM chains. However, the significant orbital magnetic 
moments appear and are also direction-dependent, except Fe. 
We also found that the Fe system has an out-of-plane mag- 
netic anisotropy whilst the Co and Ni systems have an in-plane 
anisotropy. 
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